home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / TQL1.z / TQL1
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. ____TTTTQQQQLLLL1111((((3333FFFF))))                                                            ____TTTTQQQQLLLL1111((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      TQL1, STQL1    -  EISPACK routine.  This subroutine finds the eigenvalues
  10.      of a SYMMETRIC TRIDIAGONAL matrix by the QL method.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  ttttqqqqllll1111((((nnnn,,,, dddd,,,, eeee,,,, iiiieeeerrrrrrrr))))
  15.              iiiinnnntttteeeeggggeeeerrrr          nnnn,,,, iiiieeeerrrrrrrr
  16.              ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee((((nnnn))))
  17.  
  18.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssttttqqqqllll1111((((nnnn,,,, dddd,,,, eeee,,,, iiiieeeerrrrrrrr))))
  19.              iiiinnnntttteeeeggggeeeerrrr          nnnn,,,, iiiieeeerrrrrrrr
  20.              rrrreeeeaaaallll             dddd((((nnnn)))),,,, eeee((((nnnn))))
  21.  
  22.  
  23.  
  24. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  25.      On Input
  26.  
  27.      NNNN is the order of the matrix.
  28.  
  29.      DDDD contains the diagonal elements of the input matrix.
  30.  
  31.      EEEE contains the subdiagonal elements of the input matrix in its last N-1
  32.      positions.  E(1) is arbitrary.  On Output
  33.  
  34.      DDDD contains the eigenvalues in ascending order.  If an error exit is made,
  35.      the eigenvalues are correct and ordered for indices 1,2,...IERR-1, but
  36.      may not be the smallest eigenvalues.
  37.  
  38.      EEEE has been destroyed.
  39.  
  40.      IIIIEEEERRRRRRRR is set to Zero       for normal return, J          if the J-th
  41.      eigenvalue has not been
  42.         determined after 30 iterations.  Calls PYTHAG(A,B) for sqrt(A**2 +
  43.      B**2).  Questions and comments should be directed to B. S. Garbow,
  44.      APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.